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CN ■ We discuss the accuracy of the estimation of the n identical unknown actions of 

SU(2) by using entanglement. This problem has a similar structure with the phase 
' estimation problem, which was discussed by Buzek, Derka, and Massarfl]. The es- 

timation error asymptotically decreases to zero with an order of -\ at least. 
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1 Introduction 



1.1 Estimation of unitary action 



Quantum information processing architecture is constructed with a combina- 
tion of various quantum devices. In other words, each quantum device con- 
tributes towards obtaining the output quantum state from the input quantum 
state. Hence, it is necessary to identify the quantum operation acting on each 
device. This problem is referred to as the identification of the quantum opera- 
tion; this problem usually entails the estimation of the operation based on the 
pair of the input (initial) state and the data of the measurement for the output 
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state. Generally, such a quantum operation is described by a trace-preserving 
completely positive (TP-CP) map. In this paper, we consider a noiseless case 
in which no noise appears in the quantum operation. In this case, the quantum 
operation is described by a unitary matrix. The adiabatic case is also treated 
in a similar manner. 

When the reference system is unitarily equivalent to the input system and 
the initial state is a maximally entangled state in a joint system between 
the original input system and the reference system, the final state in the 
joint system is the maximally entangled state described by the unitary matrix 
acting on the input system. Hence, by repeating this operation n times, we 
can estimate this unitary matrix by performing appropriate measurement for 
the total joint system between output system and the reference system. In the 
two-dimensional case, Fujiwara [2] showed that a maximally entangled state 
is the optimal initial state in the composite system between the single input 
system and its reference system. Ballester [3] discussed this problem in the d- 
dimensional case. When their method is applied to the case of n applications 
of this unknown operation, the estimation error decreases to zero only with 
order ^; this is the case even if the output measurement is optimized. This 
order is obtained from the accuracy of state estimation by using an identical 
state preparation. Therefore, it is worthwhile to consider the possibility of 
further improvement. 

1.2 One-parameter case: Phase estimation 

The phase estimation problem is closely related to the problem discussed in 
this paper. In the phase estimation problem, we estimate the eigenvalue e td of 
the unknown unitary matrix acting on a two-level system when we know its 
eigenvectors. The phase estimation with the fixed input state was discussed 
by Helstrom[4] first. Holevo [5] extended Helstrom's result to a more general 
framework of group covariance. The optimization problem of the input state 
for phase estimation was discussed by Buzek et al. [1] in an asymptotic setting. 
They proved that the error decreases to zero with speed ^ when we choose 
the optimal input state and optimal measurement. However, the error goes to 
zero with the order - in the usual statistical parameter estimation. Therefore, 
this unexpected result indicates the importance of the entangled input state. 

1.3 Three-parameter case: Our result 

In this study, we discuss whether such a phenomenon occurs in the estimation 
of SU(2) unitary action. In this paper, we adopt the error function d(U, U) d = 
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1 — I Tr - ^ U | 2 between the true SU(2) matrix U and the estimated matrix 
U. Then, we obtain an unexpected relation between our problem and that 
discussed by Buzek et al.[l\. Thanks to this relation, we can show that the error 
decreases to zero in proportion to ^ at least. Further, the coefficient is equal to 
7r 2 if the estimator is constructed based on this relation. It is shown that this 
bound can be asymptotically attained with no use of the reference system. 
Instead of use of the reference system, we regard a part of the composite 
system of n input systems as the tensor product of the system of interest and 
the reference system. In other words, an effect of "self-entanglement" is used 
in this method. 



2 Phase estimation: Estimation of eigenvalues with the knowledge 
of eigenvectors 



In order to demonstrate the interesting relation between phase estimation and 
estimation of SU(2) action, we briefly summarize the fact known with regard 
to phase estimation. When the eigenvectors of the unknown SU(2) matrix U 
are known, the estimation problem can be reduced to the estimation of the 
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unknown parameter 9 G [0, 2n) of the family 
because the action of U is equivalent to that of cU. 



In our framework, we assume that the tensor product matrix Uf n acts in the 
tensor product space H® n and that we can select an arbitrary initial state in 
the tensor product space Ti® n . That is, we can select the input state x G Ti® n 
and the estimating POVM M n (d9) on Ti® n with outcomes in [0, 2ir). When we 
evaluate the error between the true parameter 9 and the estimated parameter 
9 by the function sin 2 the mean error is given by 



V%(M n ,x) = [ 2 \m 2 ^— -(x\{Uf n )^M n {d8)Uf n \x) 
Jo 2 



Thus, when we fix the input state x G 7Y® n , our estimation problem can be re- 
duced to the estimation with the state family { Uf n \x)(x\{Uf n )A 9 G [0, 2?r)}. 

In order to treat this problem in a simple manner, we focus on the phase 
estimation problem in the <i-dimensional space TV spanned by the orthonor- 
mal basis ui,...,Ud as follows. We select a vector x = (xk)f =1 = Y,t=i xku k 
satisfying 



3 



EW = 1 - (!) 



fe=i 



Suppose that the state to be estimated has the form = f Uq\x){x\Uq with 

the unknown parameter 9 e [0, 27r) and the unitary matrix = f Sfe=i e* fe6, |^fc) (itfc| 
In this case, the error of the estimator M( e£0) can be be expressed as 



V e (M, x) = / sin 2 — - Tr M(d9)p e , x . (2) 

In order to ensure accuracy, it is appropriate to focus on the worst case. In 
other words, we minimize the worst error maxg T>$(M, x). This problem is 
refereed to as minimax. 

On the other hand, since the state pe, x has symmetry 



pe+e>,x — Uff/po^Ug,, (3) 
it is natural to treat the measurement M{9)d6 with the same symmetry: 

M(0 + 9') = U §l M(9)Ul. (4) 

A measurement M satisfying (4) is referred to as the covariant measurement 
and it has the form 

M T (d9)^U § TUl^, (5) 
where the Hermitian matrix T = J2k,i tk,i\ u k)( u i\ satisfies 

T > and t k , k = 1. (6) 

Holevo proved that the minimum error in the minimax criteria is equal to the 
minimum error among covariant measurements, i.e., 



min max"De(M, x) = min T>g(M,x). 

M 9 Micovariant 

This relation is referred as the quantum Hunt-Stein lemma, and was proved 
in a more general covariant setting[5]. 

By using elementary formulas of trigonometric functions, the equation 
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Jo Sm ^ u 9\uk'){w\UoU § \ui){u k \Ul — 

= fik,k>fii,i>(-$k,i — — <^fc,i— i — -£fc-i,0 (7) 

can be verified. By using this relation, we have 

Vg(M T , X) = \Xk\ tk,k - tX! (XfcX fc+ lifc+l,fc + X fc+ lX fc tfc 5 fc + l) 

2 fc=l 4 fc=l 

j d -y d-1 -y I d-1 \ 

>9 S l Xfc ! 2 ~lE (fefc+il + I^IXfel) = - I 1 - ^ |x fc ||x fe+ i| , (8) 

z fc=i 4 fc=i ^ V fc=i / 

where the inequality follows from |(it fc |T|u fc+1 )| < \J (iik\T\uk) (uk+i\T\uk+i) = 
1. Hence, the equality holds iff 



ifc,fe+i = I — n 1" ■ (9) 

\ x k\\ x k+l\ 

Since the matrix T x = f [i^jj^j l -u fc)( M d satisfies the conditions (6) and (9), 
we obtain 



umi V e (M,x) = Ul-J2\ x k\\xk+i\) ■ (10) 



Micovariant 



fc=l 



Now, we consider the estimation of the unknown unitary Ug with its multi- 
ple actions. Since the unitary matrix Uf n has eigenvalues 1, e 10 , . . . , e me , the 
application of the relation (10) to the d — n + 1 case yields 



VI* = min minmax^(M n ,x) 



(°fc))U I 2 



fe=i 



E°! = 1 . fl *>o . 

fc=0 J 



In this problem, the eigenspace Tik is not one-dimensional for the eigenvalue 
e lke . Hence, any initial state can be written as J2k x k^ k , where e k is a nor- 
malized vector in H^. Therefore, the minimum error can be expressed in the 
above form. Buzek, et al. [1] proved that the minimum error is attained by the 
coefficient a& = sin ^k+y 2 ) anc i j s almost equal to i.e., 
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In this setting, the irreducible representation space of the action of the one- 
dimensional circle S 1 = [0,2n) is the one-dimensional space spanned by u^. In 
the above discussion, we use n+1 different irreducible representation spaces. 



3 Estimation of the unknown SU(2) action 

Next, we proceed to the estimation of the unknown SU(2) action g e SU(2) 
in the two-dimensional space H. In this problem, the unitary matrix g® n acts 
in the tensor product space H,® n , and a suitable initial state in 7i® n can be 
selected for this estimation. On the other hand, Fujiwara[2] focused on the 
estimation problem where the initial state is entangled with the reference 
system Hr, in which the unknown SU(2) action g does not act. He proved 
that this method reduces the estimation error. However, he did not treat 
the entanglement between the n-tensor product space TL® n and its reference 
system. 

Now, we consider the amount of the estimation error can be reduced by use 
of entanglement among tensor product space. In other words, our framework 
has a wider choice. In our problem, we can select an initial state entangled 
between 7i® n and the reference space H% n that is unitarily equivalent to the 
original input system H® n . Hence, when we select the initial state x on the 
tensor product space H® n <S> 'H'r 1 and the measurement (POVM) M n (dg) in 
7i® n ® Hf n with the outcome in SU(2), the error is evaluated as 

V g (M, x) d ^ / d(g,g){x\(g® n ® I)*M(dg)(g® n ® I)\x). 

Here, we focus on the SU(2) action on the tensor product space H® n . Its 
irreducible decomposition is given as follows: 

d 

© ^2k+l <8> ^2(i,2fc+l 
d 

-l,2fe, 

k=l 

where Hk is the A;-dimensional irreducible space of the SU(2) action, and H n ^ 
is the corresponding irreducible space of the action of the permutation group, 
where SU(2) does not act. Note that the dimension of H. n ,k is equal to the 
number of representation spaces equivalent to Hk in the tensor product space 
U® n . Hereafter, we denote the SU(2) action on H k by V g k . 
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For simplicity, first we focus on the estimation of the SU(2) action in a sin- 
gle irreducible space Tij. Now, we select the initial state as the maximally 
entangled state x E between TCj and the reference space H.j,R that is unitar- 
ily equivalent with the space Hj. The measurement is selected as the POVM 

Mj )E (dg) = f j 2 Vj <g> I\x E )(x E \(V^ <S> iY fJ-(dg) , where /i is the invariant prob- 
ability distribution on SU(2) and the integer j 2 is the normalizing factor. By 
using Schur's lemma, we can easily verify that the total integral is constant 
times the identity matrix because the state x E is maximally entangled. Since 
(TrV^)* = Tr(yiy, the average error is calculated as 



^ w d(s,9M 2 <4l(^) t ^l4)<4l(^) t v/i4)/'W) 

jL> d{s ' S) 1 h l «-» |2 " (<i§) = L, d{1 ' §) 1 ^ ^' |2 '«» 



5 i > 2. 



(12) 



The final equation follows from the elementary calculus of trigonometric func- 
tions. This calculation seems to indicate that if we use only one irreducible 
space Tij and even if its dimension is large, the estimation error cannot be 
reduced 2 . Hence, in order to reduce the estimation error, it may be needed 
to use several irreducible spaces. 



Hereafter, we consider the case of n — 2d— 1; however, the following discussion 
can be applied to the even case. In the odd case, d distinct irreducible spaces 
exist. Hence, it is essential to use the correlation between them. We investigate 

L R 



the following subspace of H^ 2d ^ <g> Hf (2d_1) - 



m 

k=l 

Further, we denote the SU(2) representation on this space by U^ d ~ 1 . 

As demonstrated later, this problem can be treated parallel to the phase es- 
timation in which the base Uj corresponds to the maximally entangled state 
x 2 E or the vector (2k)x 2 E . We select a vector x d = (xk)(, =1 that satisfies the 
condition (1), and let the initial state be x^T 1 = f 0^ =1 XkX 2 E . In a manner 
similar to (5), the measurement is selected as 



2 This fact has been proved by Chiribella et al. [16] after the submission of the 
preliminary version of this paper. That is, they proved the optimality of this POVM 
in a more general framework. 
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Mj d _ 1 (^) d ^ f ^f- 1 E^I(2A;)4 fc )((2/)x||(^f- 1 )V(^) 
k,i 



based on a Hermitian matrix T = (t k j) satisfying condition (6). Using Schur's 
lemma and condition (6), we can verify that its total integral is an identity 
matrix. In a manner similar to (7), the equations 

^) IV C^ 4 *- 1 |a^> <a^' | (C^*- 1 ) t C^- 1 1 (2fc)a^> <(2Z)a^| (C^™- 1 ) V(d^) 

SU(2) 

=<W / su 2 %,£)TrVf (TrVf = <W'V Q<5 M - - ^ fc -i,i 

SU(2) (13) 

holds, where the final equation is derived in section 4 based on elementary 
formulas of trigonometric functions. By using this relation in a manner similar 
to (8), we have 

/ su(2) <*9,») T^r'i4f '><4r i^r 1 )*^-^) 

Y d 1 ^ 

= o X tk, k ~ 7 X ( x k-l%ktk,k-l + ^fc£fc-l4-l,fc) 
2 fc=l 4 fc=l 

a! j d_1 1 / d_1 \ 

>9 X l x fcl 2 - tE (l^k+il + l^fc+T^l) = - i - X kfclk*i+il • ( 14 ) 

2 fe=i 4 fc=i 2 V fc=i / 

The equality holds when matrix T = satisfies (9). Thus, the optimal 

error of this estimation method coincides with T)"^ 1 . That is, our problem is 
reduced to phase estimation problem by Buzek et al. When we select a suitable 
initial state and measurement in the case of large n, the estimation error is 
equal to ^ asymptotically because V d ~ t l ^ ^ ^£iT? = S- 

In the case n = 2d, the initial state is expressed as x\f == (&k=o a kxji +1 , where 

d 

the vector a' d = (a*;)fc=o has no negative element. When we select a measure- 
ment similar to the one mentioned above, the estimation error is calculated 
to be | (l — J2k=o a k a k+ij + |a . Hence, the minimum error V d pt ' satisfies 

T^tpt < T^tpt < Dipt '■ The same conclusion is obtained in the odd case. 

Next, we investigate the reference space in the odd-dimensional case. When 
we use the above method, the dimension of the reference space Hn^d-i is 2d. 
If the dimension of Ti.2d-i,2k is greater than that of Ti,2k, we can use the space 
Ti-2d-i,2k as the reference space. For k = d, the dimension of H.2d-\,2k is 1, i.e., 
is smaller than that of H,2k- However, for k < d, the dimension of H.2d-i,2k is 

( d — k ) ~~ _ k — l)' ^' e "> ^ * s S rea ^ er than that of H.2k- Hence, by replacing 
the reference space by the space H.2d-\,2kt we can reduce the estimation error to 
Dt 2 without using the reference system. Since V, 



d-2 
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an estimation error of ^ can be achieved without the reference space. This 
discussion can also be applied to the case when n is an even number. 



4 Technical details 



In the following, we investigate the derivations of equations (12) and (13). 
Since TrV~ and d(I,g) depend only on the eigenvalues of g e ie / 2 ,e~ l6 / 2 , the 
invariant distribution /j, has the form fi(dg) = - sin 2 | d9 sin (fiid(j)id(f)2 in the 
following parameterization: 



9 = 



cos 0i sin (pie 1 ^ 2 
^ — sin (j)ie^ 2 cos 0i 

Hence, the relation 



/ / P^/ 2 n W cos 0i sin 0i 



e™'- 
e'^ 2 



y — sin 0i e 92 cos 0i J 



f f(6Mdg) = r M sin 2 9 de 
Jsv<2) Jo n 2 



/SU(2) 



(15) 



holds. Because d(I, g) = sin 2 § and Tr V? 
(15), we have 



= Y?i=i^ 1 3+2 ls>e an d by applying 



SU(2) 



d(I,g)\TrVn 2 v(dg) 



27r -sin 4 -^^^-^)^ ^ -"-^ 

7T 



r 2w 1 /3 1 \ / J \ 



When j > 2, this integral is equal to 



2 ^ 1 /3 
4^ V2 



2 cos 6> + - cos 20^ {j + 2(j - 1) cos + 2(j - 2) cos 20) d0, 
which implies (12). When j = 1, we have 



r2n 1/3 1 \ 1 

/ — (--2cos0+-cos20 )<$=-, 
io 4tt V2 2 / 2' 



which implies (12) 
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With regard to (13), we can calculate 

In the above integral, the term (j2 2 =i e^' -2 ^) ) (j2 2 =i e 1 ^ -2 ^) 6 ^ can be ex- 
panded into many terms; however, all other terms except the constant term, 
2cos6>, and 2 cos 2$ vanish. That is, the coefficients of only the above three 
terms are important. 

In the following, we calculate this integral in the three cases: (i) k = k', (ii) 
| k — k' | > 1, and (iii) \k — k'\ = 1. Case (i) has already been calculated 
in (12). In case (ii), these three coefficients coincide. Thus, the integral is 
equal to 0. Next, we proceed to the case (iii). When k! = k + 1, the term 

(E?=i e i{ - l -^ d ) (E/='i e^" 2 ^) can be expanded to A; + 2k cos 9 + 2(k - 
1) cos 29 + • • • . Thus, the integral is equal to — |. Similarly, we can show that 
the integral is equal to — \ when k' — k — 1. Therefore, we obtain (13). 



5 Concluding remark 



In this paper, we derived a remarkable relation between the estimation of the 
SU(2) action and phase estimation. By using this relation, we showed that 
the optimal estimation error is less than ^. The essence of this relation lies 
in the relation between the two similar relations (7) and (13). Indeed, since 
sin 2 9 = 1 — | e t8 / 2 +e 10 2 12^ ^ e cog ^ f unc tion of phase estimation is equal to 

1 — \X2(9 — 9)\ 2 , where X2(9) is the half of the character of the two-dimensional 

. Hence, the essence of (7) is the equation 



representation e 



^ e'^ 2 J 



L 



2vr 

2 



(1 - \X2(9)\ 2 ) X k (0)x l (0)d9 = -8 Kl - - -5 k ^, (16) 

where rf(0) is the character of the one-dimensional representation of e %e i— > e de . 

On the other hand, when we let x l {d) assume the character of the /-dimensional 
irreducible representation of SU(2) over dimension I, the essence of (13) is 



/su ( 2) (1 " \^(9)\ 2 )x 2k (9)x 2l (9Mdg) = \s kJL - - (17) 
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Therefore, the main reason for the relation between the two estimation prob- 
lems is the character formulas (17) and (16). In other words, the integrals with 
regard to characters are the essence of phase estimation as well as that of the 
estimation of SU(2) action. Physically, the square-error estimators in both 
scenarios can be realized by the effective use of interference between different 
irreducible representations. 

In this paper, when we consider the irreducible space 7i n> k of the action of 
the permutation group as the reference space of the irreducible space Tik 
of SU(2) action, we proved the existence of a square-error estimator in the 
estimation of the SU(2) action without use of the real reference system. Since 
the irreducible space 1~in,k of the action of the permutation group has a greater 
dimension than the corresponding irreducible space Hk of the SU(2) action, 
except for the exceptional cases, the correspondence between the irreducible 
space TL n ^k and the reference space of the corresponding irreducible space TCk 
is not unique. Hence, it is desirable to obtain a more physically realizable 
correspondence [13]. This is an interesting future study. 

Finally, we should remark on the relation between our group covariance ap- 
proach and the Cramer-Rao approach. In the latter approach, we focus on the 
Cramer- Rao-type bound, i.e., the minimum weighted sum of mean square er- 
rors under locally unbiased conditions. As discussed in Fujiwara and Imai[7], 
we often seek an input state that minimizes this bound. Even though this 
bound decrease to zero with an order of \, we cannot conclude that there 
exists a sequence of estimators whose weighted sum of minimum mean square 
errors decreases to zero with an order of A. In the state estimation, there exists 
a sequence of estimators whose weighted sum of minimum mean square errors 
reduces to zero at rate — , where C is the Cramer- Rao- type bound. This is be- 
cause an adaptive estimator attaining bound C can be selected as follows[8,9]. 
In the state estimation, we can select a set of neighborhoods in which the 
minimum weighted sum of mean square errors can be approximated to the 
Cramer-Rao type bounds within bounded differences. Hence, if we choose the 
first estimator whose estimate belongs to the above neighborhood with an 
exponentially small error probability, the adaptive estimator satisfies the re- 
quired condition. However, the situation of the estimation of the SU(2) action 
differs from that of the estimation of state. In the estimation of the SU(2) 
action, such neighborhoods depend on the number n because the state fam- 
ily depends on the input state and the required neighborhoods depend on 
the state family. In other words, there is a possibility that the radius of the 
neighborhoods reduces to zero in proportion with the number n of actions. 
Hence, we cannot directly obtain the estimation error based on the optimal 
Cramer-Rao-type bound. 

On completion of this research, the author found that the same results were 
obtained by two other groups[10,ll,12]. However, this approach is different 
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from those employed by them because it is based on the notable relation 
between the SU(2) estimation and the phase estimation. Further, Chiribella 
et al.[lQ] proved the optimality of the error | ( 1 — J2t=i \ x k\ l^fe+il) while this 
is not proved in this paper. They also elucidated out that character integrals 
and interference of irreducible representations are very useful for finding the 
optimal estimation of an unknown group transformation not only for U(l) and 
SU(2) but also for arbitrary groups. 

Further, on completion of this paper, the author also found Acin et aVs paper 
[15]. They almost mentioned that the optimal initial state has the form x^ -1 . 
However, they did not derive the equation corresponding to (14). In addition, 
Rudolph and Grover [14] discussed a similar problem from a computational 
viewpoint. 
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